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ABSTRACT 

The effect of stratospheric aerosols on the earth's monthly zonal radiation 
balance is investigated using a model layer consisting of 75% H 2 SO^, which is the 
primary constituent of the background aerosol layer. The reduction in solar 
energy absorbed by the earth-atmosphere system is determined through the 
albedo sensitivity, defined here as the change in albedo per unit mid-visible 
optical depth of the aerosol layer. The optically thin approximation is used in 
conjunction with the Henyey-Greenstein phase function for scattering to simplify 
computations. Satellite derived planetary albedos are used as the frame of 
reference about which the change in albedo is computed. An infrared radiative 
transfer model is used to estimate the increased greenhouse effect attributed to 
the aerosol layer. The infrared heating tends to compensate for the albedo 
effect in altering the radiation balance. The results indicate that the dominant 
influence of the thin model stratospheric aerosol layer is an increased reflection 
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of solar energy all over the globe except for the polar-winter region, but the 
change in the radiation balance is seen to be imiform and small equatorwards of 
The largest perturbations occur in the spring and fall in tlie polar regions 
with the equator-to-pole radiation balance increasing by 6, 5-7.0 W/m- for an 
aerosol mid-visible optical depth of 0.1. 
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ALBEDO ENHANCEMENT AND PERTURBATION 
OF RADIATION BALANCE DUE TO 
STRATOSPHERIC AEROSOLS 

1. INTRODUCTION 

The role of stratospheric aerosols In enhancing the reflection of solar ra- 
diation has been investigated recently by numerous authors (Harshvardhan and 
Cess, 1976; Pollack, et al. , 1976; Cadle and Grams, 1975; Pinnick, et al. , 

1976; Herman, et al. , 1976). Model calculations by the above and also by 
Luther (1976) have shown that the dominant climatic effect of tliese aerosols is 
enhanced solar reflection with increased infrared opacity playixig a smaller role. 
However, Herman, et al. , (1976) and Luther (1976) have also pointed out that 
the degree of albedo enhancement is a strong fmiction of the albedo of the imder- 
lylng surface. Broadly spealdng, tliis imdeiiying surface may be divided into 
two classes, viz. cloud-free portions of the earth-atmosphere system and the 
cloudy portion. Tins classification does not imply a Ihgh and low albedo differ- 
entiation because the effective albedo of a cloud layer above a snow or ice cov- 
ered surface may be lower than that of the surface alone (Dave and Braslau, 
1975). However, with the above exception, the effective albedo of the cloud- 
free regions is between 0. 15 and 0.3 whereas that of cloudy regions is greater 
than 0. 4. In this discussion, the reflection and absorption in tlie atmosphere 
is included in computing the effective albedo, an implicit assumption being that 
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the reflectivity of the atmosphere above the stratospheric aerosol layer may be 
neglected. 

The purpose of the present study is to examine the albedo sensitivity of the 
earth-atmosphere system to a stratospheric aerosol layer under different as- 
sumptions regarding surface and cloud albedos and solar zenith angle averaging 
procedures. Work to date encompasses a broad range of simplifying assump- 
tions. Harshvardhan and Cess (1976) considered an albedo of 0. 3 for the unper- 
turbed earth- atmosphere system. Pollack, et al. , (1976) assumed a ground 
albedo of 0. 1 and 50% cloud cover with a detailed radiative transfer scheme to 
yield aerosol layer perturbations, while Herman, et al. , (1976) assigned ground 
albedos to 10° latitude belts in tlie Nortliern Hemisphere with seasonal variations 
and assumed a 50% cloud cover of reflectivity equal to 0. 5. The marked changes 
in albedo sensitivity obtained by Herman, et al. , with changing seasons at high 
latitudes, due to changes in the mean solar zenith angle and surface albedo, 
suggests that incorrect inferences may be drawn if global average surface 
albedos and cloud cover are assumed in computing albedo enhancement. 

In this work the albedo of the cloud-free portion of the earth-atmosphere 
system is presented in 10° latitude belts for each month of the year following the 
recommendation of Budyko (1974). An empirical relation for the effective cloud 
reflectivity is derived from calculations of Dave and Braslau (1975). With tliese 
two scjts of albedos as a base, the global albedo change in the presence of a 
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model stratospheric aerosol layer is computed for climatological cloud cover 
fractions. The pi'ocedure is extended to include the zonal variation of albedo 
sensitivity by month, based on satellite derived planetaj.’y albedos. With the 
inclusion of the effect of the stratospheric aerosol layer on the infrared radia- 
tion leaving the top of tiie atmosphere, the net energy deficit on a zonal basis is 
obtained. This quantity may be used as an input to energy-balance climate 
models used for investigating the climatic effects of stratospheric aerosols. 

2. ALBEDO SENSITIVITY 

The system albedo, cvp, is commonly expressed as 

ttp = ttcC + as(l -C) (1) 

where ag is the cloud-free or clear sl^y albedo and tx^ the albedo with cloud 
cover; C is tlie fractional cloud cover. The expression for ap may be used on 
a global basis or for individual latitude belts (Budyko, 1974; Cess, 1976). 

For the remainder of tlie text «p will be used to denote the albedo for each 10° 
latitude belt, witli the hemispheric planetary albedo, A, computed as 

A = - / ctp(^)QW)cos^d0 (2) 

^ •'0 

where 0 is the latitude, Q the latitudinal distribution of incoming solar radiation 
and S the solar constant. The integration of Equation (2) is replaced by a 
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summation over the nine hemispheric latitude belts. The perturbation to this 
planetary albedo with the introductioii of a stratospheric aerosol layer enhance- 
ment (as following a major volcanic eruption) is related to the optical depth 
of the aerosol layer. The albedo sensitivity may be defined as dA/dr where the 
optical depth, r, is computed at some reference wavelength, usually the mid- 
visible wavelength of 0. 55 pm. Since the stratospheric aerosol layer is optically 
thin (r«i), the change in albedo is directly proportional to the change in op- 
tical depth and the albedo sensitivity may be written as AA/Ar (Russell and 
Hake, 1977), In succeeding sections, models of the planetary albedo and reflec- 
tivity of the aerosol layer are postulated and the albedo sensitivity computed. 

3. SURFACE AND CLEAR SKY ALBEDO 

The fraction of solar radiation incident on a surface that is reflected by it 
is the surface albedo Using the mean value of albedo for various natural 

surfaces, Budyko (1974) has suggested appropriate values of land albedos to be 
used for different climatic and seasonal regimes. The values vary from 0. 8 for 
stable snow cover in high latitudes to 0, 13 for steppe and forests. The albedo of 
a water surface to be used for the ocean fraction of latitude belts has also been 
given by Budyko based on the work of Sivkov (in Kondratyev, 1975), As there is 
a strong zenith angle dependence of the direct beam albedo, there is a marked 
seasonal variation at naid- and high-latitudes. The albedo of both land and ocean 
change dramatically with snow and ice cover and this seasonal variation must 
also be incorporated in the computation of global albedo. In the present study. 
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data on snow and ice coyer compiled bj’^ Curran, et al. (1978) were used to 
generate the seasonal march of albedo at liigh latitudes. The surface albedos 
for land and ocean computed after identifying the surface type and incorporating 
snow and ice cover are presented in Tables 1 and 2, respectively. 


The clear-sky albedo, , is the effective albedo of a cloud-free atmosphere 
over the land and ocean surface identified by the albedos, q-sq , given in Tables 
1 and 2. It is possible to mal^e a detailed radiative transfer calculation averaged 
over the solar spectrum to obtain ttg ; however, for our purposes, a simple para- 
meterization that takes into account molecular and particle scattering as well as 
gaseous absorption in the atmosphere is sufficient. One such relation may be 
obtained from the model of Braslau and Dave (1973) for midlatitude summer 



water vapor and ozone distributions and an absorbing troposp]i|^*-lcj aerosol layer. 
Choosing their Model Cl, the planar reflectivity of a-€iear sky over a black sur- 
face may be vTitten as 

as(R = 0,G) = 0.05 + 2.4 x 10-^02 + 3.2 x (3) 


where 0 is the solar zenith angle in degrees and R is the reflectivity of the under- 
lying surface. To extend this relation to reflecting surfaces, a multiple reflection 
model is used which yields (Coaldey and Chylek, 1975; Wiscombe, 1975) 

RT Trl 

a-,(R,0) = a,(R = O,0) + (4) 

1 - Raj(R= 0) 
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Month 

Latitude '^.^ 

B 

F 

N 85 

— 

— 

75 

— 


65 

— 

0.80 

55 

0.70 

0.70 

45 

0.70 

0.70 

35 

0.20 

0.20 

25 

0.25 

0,25 

15 

0.20 

0.20 

N 5 

0.10 

0.10 

S 5 

0.10 

0.10 

15 

0.20 

0.20 

25 

0.20 

0.20 

35 

0.20 

0.20 

45 



55 




65 

75 

85 


S 


0.20 

0.80 

0.80 


0.20 

0.80 

0.80 


M 


A 


S 


O 


N 









No 

Land 

M 

M 


— 

— 

0.80 

0.50 

MSk 

mm 

0.50 

0.80 

— 

0.30 

0.20 

19 

19 

0.15 

0.30 

0.80 

0.15 

0.15 

iig 

Hg 

0.15 

0.30 

0.70 

0.15 

0.15 

0.15 

0.15 

0.15 

0.20 

0.40 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

0.10 

0.10 

0.10 

0.10 

0.10 

0.10 

0.10 

0.10 


19 

0.10 

0.10 


WKM 

0.20 



0.20 

0.20 

91 

9^9 

0.20 

0^9 

0.20 

0.20 

0.20 

mm 

9^9 

0.20 

0.20 

0.20 

0.20 

0.20 

0.20 

9^ 

No 

Land 






0.80 

— 

— 

0.80 

0.80 

0.40 

0.20 

— 

— 

— 

— 

0.80 

0.80 

0.80 

„ ^ 

— 

— 

— 

0.80 

0.80 

0.80 





















Table 2 


Albedo of Ocean and Sea Ice 


Month 

J 

F 

' M 

A 

M 

J 

J 

A 

S 

0 

N 

D 

N 85 

— 

— 

0.80 

0.80 

0.80 

0.70 

0.70 

0.70 

0.70 






75 


— 

0.80 

0.80 

0.80 

0.70 

0.70 

0.70 

0.70 

0.80 


— 

65 

~ 

0.31 

0.28 

0.26 

0.25 

0.22 

0.19 

0.17 

0.17 

0.20 

0.26 

— 

55 

0.21 

0.18 

0.15 

0.13 

0.11 

0.09 

0.08 

0.08 

0.09 

0.12 

0.17 

0.20 

45 

0.13 

0.10 

0.09 

0.07 

0.07 

0 . 06 

0.06 

0.07 

0.08 

0.09 

0.12 

0.14 

35 

0.10 

0.09 

0.08 

0.07 

0.06 

0.06 

0.06 

. 0.06 

0.06 

0.07 

0.09 

0.10 

25 

0.08 

0.07 

0 . 07 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.07 

15 

0,07 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.07 

N 5 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

S 5 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

0.06 

15 

0,06 

0.06 

0.06 

0.06 

0.06 

0.07 

0.07 

0.06 

0.06 

0.06 

0.06 

0.06 

25 

0.06 

0.06 

0.06 

0.06 

0.07 

0.08 

0.08 

0 . 07 

0 . 07 

0.06 

0.06 

0.06 

35 

0.06 

0.06 

0.06 

0 . 07 

0 . 

0. 10 

0.10 

0.09 

0.08 

0.07 

0.06 

0.06 

45 

0.06 

0.07 

0.08 

0.09 

0.12 

6. 14 

0.13 

0.10 

0.09 

0.07 

0 . 07 

0.06 

55 

0.08 

0.08 

0.09 

0.12 

0.16 

0.18 

0.19 

0.16 

0.13 

0.12 

0.10 

0.07 

65 

0.16 

0.12 

0.15 

0.22 

0.30 

— 

— 

0.36 

0.35 

0.34 

0.30 

0.21 

75 

S 85 

0.70 

0.70 

0.70 

0.80 

No 

Water 



0.80 

0.80 

0.80 

0.70 





























where Tc ajad Tj are the planar and diffuse transmission respectively and 
(R = 0) is the diffuse or global albedo over a black surface. The atmospheric 
absorptivity used to compute T^ and Tj has been assumed to be 15% which is an 
average value based on Braslau and Dave’s computations. The variation of ab- 
sorption with solar zenith angle is not very severe and has been ignored. The 
results of tliis model are compared with Braslau and Dave (1973) in Fignre 1. 

The correspondence is excellent for both low and high values of R; the differ- 
ences are least at the higher zenith angles. The mean solar zenith angle, even 
at the equator is always greater than 50°. Hence at high latitudes, where R is 
large due to snow and ice cover, 6 is much larger so the parameterization is 
adequate. 

Using the above parameterization, one can generate the clear-sky albedos 
from surface albedo values and these are given in Table 3. The final column in 
the table is the insolation weighted mean clear-sky albedo of each latitude belt. 
This is the value of that should be used in mean annual global energy balance 
climate models as it is a measure of the total solar radiation reflected back to 
space by the clear sky portion of each latitude belt on an annual basis. The 
values of clear sky albedo are compared with the minimum albedo obtained from 
29 months of satellite observations by Vonder Haar and Ellis (1975) in Figure 2. 
Although sampling by various satellites occurred at just one local time during 
daylight hours, the 29-month set includes late morning, near noon and early 
afternoon orbits; so the results may be considered to correspond to a mean solar 
zenith angle. 
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CLEAR SKY ALBEDO 
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MULTIPLE REFLECTION 
MODEL 



0 L. I \ I I \ I I t I 

0 10 20 30 40 50 60 70 80 90 

SOLAR ZENITH ANGLE deg. 

Figure 1. Clear-Sky Albedo Using Multiple Reflection 
Model Compared mth Braslaii and Dave (1973) for 
Surface Reflectivity of 0, 0.1, 0. 3, 0. 6 and 0.8 
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Table 3 


Cloud- Free Albedos 


Latitude 

Month 

■ 

F 

M 

A 

M 

J 

J 

A 

S 

0 

N 

D 

Insolation 

Weighted 

Mean 

N 

85 

— 



0.64 

0.63 

0.61 

0.54 

0.54 

0 . 56 

0 . 58 







0 . 58 


75 

— 

— 

0.63 

0.62 

0.61 

0.50 

0.49 

0.50 

0.54 

0.64 

— 

— 

0.54 


65 

— 

0.55 

0 . 53 

0.51 

0.28 

0.23 

0.21 

0.21 

0.23 

0.33 

0.56 

__ 

0.30 


55 

0.46 

0.43 

0.40 

0.25 

0.18 

0.17 

0.17 

0.17 

0.19 

0.28 

0.45 

0.46 

0.24 


45 

0.40 

0.37 

0.26 

0.16 

0.15 

0.15 

0.15 

0.15 

0.17 

0.21 

0.30 

0.40 

0.20 


35 

0.22 

0.20 

0.17 

0.16 

0.15 

0.15 

0.15 

0.15 

0.16 

0.18 

0.20 

0.22 

0.17 


25 

0.19 

0.18 

0.17 

0.16 

0.16 

0.16 

0.16 

0.16 

0.16 

0.17 

0.19 

0.20 

0.17 


15 

0.16 

0.14 

0.13 

0.13 

0. 13 

0.13 

0.13 

0.13 

0.13 

0.14 

0.15 

0.16 

0.13 

N 

5 

0.13 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0.12 

0.13 

0.12 

S 

5 

0.12 

0.12 

0.12 

0.12 

0.12 

0.13 

0.13 

0.12 

0.12 

0.12 

m 

m 

0.12 


15 

0.13 

0.13 

0.13 

0.14 

0.15 

0.16 

0.16 

0.14 

0.13 

0.13 


Ha 

0.13 


25 

0.13 

0.13 

0.13 

0.15 

0.17 

0.18 

0.18 

0.16 

0.15 

0.13 

Ha 

0.13 

0.14 


35 

0.13 

0.14 

0.16 

0.18 

0.19 

0.19 

0.18 

0.15 

0.14 

0.12 

0.12 

0.12 

0.14 


45 

0.13 

0.15 

0.18 

0.22 

0.24 

0.23 

0.19 

0.17 

0.14 

0.13 

0.12 

0.12 

0.15 


55 

0.15 

0.15 

0.17 

0.22 

0.28 

0.30 

0.30 

0.25 

0.21 

0. 18 

0.16 

0.14 

0.18 


65 

0.21 

0.19 

0.24 

0.34 

0.41 

— 

— 

0 . 41 

0.38 

0.33 

0.28 

0.24 

0.26 


75 

0.59 

0.60 

0.61 

0.64 

— 

— 

_ 

— 

0.63 

0.62 

0.61 

0.58 

0.60 

S 

85 

0.61 

0.62 

0.64 

— 

— 


— 

— 

0.64 

0.63 

0.61 

0.60 

0.61 
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LATITUDE 

Figure 2. Model Clear-Sky Albedos Compared with 
Minimum Albedos from Vonder Haar and Ellis (1975) 


11 . 



Viewing geometry and the presence of clouds over snow and ice make observa- 
tions near the poles unreliable; so only the portion of tlie globe between latitudes 
65°N and S should be considered. Substantial differences between the model and 
the minimum albedo exist at high latitudes in the Southern Hemisphere and near 
the Equator. The former is a region of extensive cloudiness (up to 80% average 
cloud cover; van Loon, 1972), and it is possible there were no totally cloud free 
satellite viewings recorded, resulting in the 5-10% positive discrepancy. Since 
tiie atmospheric model uses one set of radiativel 3 '- active constituents to repre- 
sent the whole globe, the parameterization is adequate. 

4. CLOUDY-SKY ALBEDO 

The spectrally averaged reflectivity of an isolated cloud layer is a function 
of the physical properties of the cloud and the zenith angle of the incident solar 
radiation. The albedo of a cloudy sky further depends on the reflectivity of the 
underlying surface. Since the number of variables involved in specifying a 
cloud completely for a radiative transfer computation is large and the nature 
of cloud cover and underlying surface is variable, parameterization of the 
reflectivity is necessary. In the present work, the cloudy-sky albedo is 
assumed to depend on the mean solar zenith angle and reflectivity of the under- 
lying surface. An empirical relation betv^een the effective albedo of a cloudy 
sky over a non-reflecting surface, given by Dave and Braslau (1975) is extended 
to yield cloudy-sky albedos using a multiple reflection model as in Wiscombe 
(1975). 
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Dave and Braslau (1975) have presented results for a detailed atmospheric 
model with midlatitude summer water vapor and ozone distributions and two 
tropospheric aerosol concentrations, and also a non-absorbing cloud layer 
between 3 and 4 km. For parameterization, their model with average aerosols 
(Cl-ST) has been chosen and an excellent fit to the planar reflectivity over a 
black surface is obtained with 

ac(R=0,6)) = 0.194 + 4.9 x 1O'502 ( 5 ^ 

where d and R have been defined previously. To extend tliis relation to reflecting 
surfaces. Equation (4) is used, again vdth 15% absorption. The results of tliis 
model are compared with Dave and Braslau (1975) in Figure 3. Again the 
parameterization is found to be adequate for the range of II and d encountered. 

The cloud model used by Dave and Braslau has an optical depth of 3. 35 at 
0. 55 pm and a diffuse albedo of 33.2%. This model (henceforth called Model C) 
has been used to calculate the cloudy-sky albedo, , for each latitude belt and 
month of die year and is presented in Table 4. The insolation weighted annual 
mean cloudy-sky albedo is shown in Figure 4 along with tlie values of oc^ deduced 
by Cess (1976) using satellite measurements of the -global albedo and mean 
annual cloud cover. The model is seen to be an imderestimate at all but equa- 
torial latitudes, indicating tlaat the optical depth used in tlie model is too low. 

Dave and Braslau (1974) have also presented planar albedos of isolated cloud 
layers of different optical depths. Choice of an optical depth of 6 and 
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Figure 3. Cloudy- Sky Albedo Using Multiple Reflection 
Model Compared with Dave and Braslau (1975) for 
Surface Reflectivity of 0, 0.1, 0. 3, 0.6 and 0.8 






Table 4 


Cloudy-Sky Albedos (Model Q 


Latitude 

Month 

B 

F 

M 

B 

M 

B 

B 

B 

s 

0 

N 

D 

Insolation 

Weighted 

Mean 

N 

85 

— 



0.72 

M 

0.68 

0.63 

0.65 

0.66 

0.69 

_ 

_ 

_ 

0. 66 


75 

— 

— 



0.69 

0.62 

0.62 

0.63 

0.65 

0.72 

— 

— 

0.65 


65 

— 

0.67 

0.64 

n 

0.48 

0.46 

0.45 

0.46 

0.49 

0. 55 

0.67 

— 

0. 51 


55 

0. bi 

0.58 

0.56 

0.46 

0.42 

0.41 

0.41 

0.42 

0.45 

oi 51 

0. 59 


0.46 


45 

0.56 

0.54 

0.47 


0.39 

0.38 

0.39 

0.39 


0.46 

0. 51 

0. 57 

0.43 


55 

0.46 

0.44 

0.41 

0.38 

0.38 

0.38 

0.38 

0.38 


0.42 

0.45 

0.47 



15 

0.43 

0.42 

EXE] 

0.38 

0.37 

0.37 

0.37 

0.38 

0.39 

0.42 

|X3 

0.44 



15 

ESI] 

0.39 

0.37 

0.36 

0.16 

ESI 


0.36 

0.37 



0.41 

0.37 

N 

5 

0.16 

0.17 

0.35 

0.35 

0.37 

0.37 

m 

0.36 

0.35 


0.17 

0.38 

0.16 

8 

5 

0.18 

0.15 

0.15 

0.16 

fm 

fm 

0.38 

0.37 

0.35 


0.18 

0.17 

0.16 


15 

0.18 

0.18 

0.18 

0.18 

EXT] 



0.18 

0.17 


0.18 

0.18 

0.17 


15 

0.18 

0.18 

0.17 

0.19 

o.u 


0.42 

0.41 

0.18 

0.17 

0.18 

0.18 



15 

0.18 

0.17 

0.18 

0.41 

0.44 

0.48 

0.45 

o.a 


0.18 

0.17 

0.18 

0.19 


U 

0.18 

0.18 

0.41 

0.45 

0.48 

0.48 

0.48 

0.48 


Mil 

0.18 

0.18 

0.41 


55 


0.41 

0.44 

o.a 

0.51 

0.54 

8.51 

0.50 


0.43 

0.41 

0.40 

0.41 


81 

0.45 

0.45 

0.48 

0.58 

0.80 

B 

— 

0.80 


0.52 

rwn 

0.47 

0.48 


fi 

0.88 

0.88 

0.88 

0.71 

B 

B 

— 

— 

0.71 

0.70 

0.88 

0.87 

0.88 

8 

88 

8.88 

0.78 

8.71 

— 

H 

H 

— 

— 

0.71 

0.70 

KU 

m 

0.88 














































CLOUDY SKY ALBEDO 



NORTH SOUTH 


LATITUDE 

Figure 4. Cloudy-Sky Albedos for Three Models 
Compared with Values Deduced by Cess (1976) 
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assumption that the cloud layer is non-absorbing, results in an equation like 
Equation (5): 

ac(R=O,0) = 0.3 + 4.25 x lQ-502 (6a) 

henceforth called Model A. Another model albedo intermediate between that of 
Equation (5) and Model A is assumed to be 

ac(R=O,0) = 0.25 + 4.625 x 10-502 (6b) 

and is called Model B. The mean annual results have been presented in Figure 
4 along with the model given by Equation (5) and Cess (1976). 

The highest reflectivity is given by Model A and is seen to be an overesti- 
mate throughout while Model B provides good correspondence at mid and high 
latitudes. The three cloud models have been used in Equation (1) to give the 
latitudinal distribution of mean global albedo using the clear-sky albedos of 
Figure 2 and cloud cover fractions, C, from London (1957) for the N. Hemi- 
sphere and van Loon (1972) as reported by Cess (1976) for the S. Hemisphere. 
The results are presented in Figure 5. Also marked are values of ctp from 
satellite measurements (Ellis and Vender Haar, 1976) and from the Nimbus-6 
ERB experiment (Jacobowitz, et al. , 1978). Except for latitudes poleward of 
65°, the variation is seen to be simulated quite well by Model B at midlatitude. 
The N. Hemisphere planetary albedo, A, defined by Equation (2) for the three 
models is 0. 323 for Model A, 0. 308 for Model B and 0. 291 for Model C. 
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LATITUDE 

Figure 5. Planetary Albedos Compared with Satellite 
Derived Results 
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It may be noted tliat the planetary albddt? derived from the satellite resiilts of 
Ellis and Vonder Haar is 0. 306. The three models are seen to cover the range 
of reasonable values of the cloudy-sky albedo and the planetary albedo and as 
such may be used as the basis for computing the albedo sensitivity. 

5. STRATOSPHERIC AEROSOL LAYER 

There is a persistent tenuous layer of submicron size aerosols in the 
stratosphere that is composed primarily of an aqueous solution of sulfuric acid. 
The number concentration of this aerosol layer increases after major volcanic 
events, and as the residence time of these aerosols is of the order of a year or 
more, this enlianced concentration persists for up to tlmee years. Recent 
measurements (Hofmann, et al. , 1975; Rosen, et al. , 1975) dux’ing a period 
free of major volcanic activity have given a picture of the worldwide distribution 
of the background stratospheric aerosol layer. For purposes of radiative trans- 
fer studies, the quantity of interest associated with the backgromid layer is the 
optical depth which is 0. 005 at 0. 55 [im, a value that may be deduced from 
Hofmaim, et al. 's results with an assumption regarding size distribution and 
refractive index of the aerosols. This illustrates the tenuous nature of the back- 
ground aerosol layer which may safely be neglected from most atmospheric 
radiative transfer computations. 

Following a volcanic eruption in which the ejected material penetrates into 
the stratosphere, there is a very significant enliancement of the aerosol layer. 
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Even after the larger particles that form part of the volcanic debris have fallen 
back into the troposphere, the aerosols formed from the gaseous ejection remain. 
Initially this enhanced concentration is limited to the vicinity of the explosion but 
in a matter of months may spread globally (Dyer and Hicks, 1968). The enhance- 
ment of the layer may be 10 or 20 fold in the latitude belt containing the explosion 
(Volz, 1970) and this can definitely be expected to have radiative effects. 


One result of this severe enhancement is the additional reflection of solar 
energy back to space, and previous work mentioned in Section 1 shows that this 
indeed is the major radiative effect. However, tlie magnitude of this effect is 
crucial in determining the radiative perturbation that may be expected from an 
enhancement. The albedo sensitivity introduced in Section 2 is a measure of 
the increased solar reflectivity. A radiative model of the aerosol is needed to 
compute the reflectivity of the aerosol layer. Mie theory for spherical particles 
may be used to compute radiative parameters. The aerosol composition is 
assumed to be 75% H2SO4 (Rosen, 1971; optical constants tabulated in Palmer 
and Williams, 1975) and a modified gamma size distribution (Shettle and Fenn, 
1976), normalized to 1 particle/cm^ is adopted: 


dn(r) 

dr 


324r exp(-l Sr) 


( 7 ) 


with dn(r)/dr the number of particles per cm^ with radii between r and r + dr, 
where r is in microns. At the mid-visible wavelength of 0. 55 pm, the radiative 
parameters are = 1,1 x 10''^ km"^, c5 = 1 and g = 0.73 with and the 
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extinction and scattering coefficients respectively, cj the single scattering al- 
bedo and g the asymmetry parameter. The optical depth of a uniform layer of 
number concentration 1 cm”^ is ~ ^ex ^ where h is the height or depth of the 
layer. A typical background distribution at mid-latitudes is given in Table 5. 


6. AEROSOL REFLECTIVITY 


As in the case of a cloud deck overlying a reflecting surface, a multiple 
reflection model may be used to determine the albedo sensitivity of the strato- 
spheric aerosol layer. For this, the reflective property of the isolated aerosol 
layer has to be determined based on the model introduced in Section 5. Here 
all computations are made at a wavelength of 0. 55 fxm.; however, radiative 
properties may be averaged over the solar spectrum. The reflectivity of the 
layer at 0. 55 /iin is roughly 25% greater than a solar averaged value 
(Harshvardhan and Cess, 1976). 

Since the optical depth of the background aerosol layer is only 0. 005 
and even after major volcanic events is normally less than 0. 1, the optically tliin 


approximation is used to compute the reflectivity of the layer. In tliis limit the 
planar albedo, a, is given by (CoaMey and Chylek, 1975) 



where r is optical depth, Mq the cosine of the angle of incidence and p(M, m' ) is 
the azimuth independent part of the scattering phase function. The global albedo, 
a, is 
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Table 5 


Background Distribution of Stratospheric Aerosols at Mid-Latitudes 
(based on Shettle and Fenn, 1976) 


Altitude, km 

Pexf lair 

1 (at 0. 55 jam) 

Spring-Summer 

Fall-Winter 

11 

8.0 X lO-'l 

6.0x10-“-™^""'^ 

location 

12 

5.5 

5.5 

13 

5.0 

5.5 

14 

^ g tropopause 

location 

5.5 

15 

4.0 

6.0 

16 

4.0 

6.5 

17 

4.5 

6.0 

18 

5.0 

5.5 

19 

5.2 

5.0 

20 

5.3 

4.0 

21 

5.0 

3.2 

22 

3.5 

2.7 

23 

2.5 

2.0 

24 

2.0 

1.7 

25 

1.6 

1.4 

Total Optical Depth 
25 

"^O-SS/Lim ~ S /^ex 

0.0054 

0.0057 








a = cor 


*'0 •'0 


dM' P(M,-m') 


( 9 ) 


It may be noted that both a( Mq ) and a are directly proportional to the optical 
depth in this limit. The scattering phase function is computed from the model 
parameters mentioned in Section 5 and the integrations in Equations (8) and (9) 
may be carried out for selected solar zenith angles following the procedure of 
Wiscombe and Grams (1976). The use of the Henyey-Greenstein phase function 
(Hansen, 1969), which is 


pH ccn - 


1 - 

(1 + g2 -2gn^/2 


( 10 ) 


where g is the asymmetry parameter and ^ the cosine of the scattering angle, 
introduces only small errors in flux and is used here because closed form 
solutions can be obtained for a(ju) and a. These are 


COT 

a(A^o) " ~ 

Mo 


^0 - — Mo - 


^22+J P2C+]fMo) 

£=] 


~ (-l)(-3)(-5) . . . {-2e+l) 

J2K + 2)(2e)(2e - 2) .... 2 


( 11 ) 




^2 2+ 1 


~ (-l)(-3)C-5)...(-2e-M) ~|~ I 

_(2fi + 2)(2e)(2e-2) 2 J j 


( 12 ) 


where c5o= 1> = 3g/4 and c52g+] = (42 4-3)g^^'^^ and Pg is the Legendre 

Polynomial of order 2 . 


The above equations give the reflectivity of an isolated aerosol layer for 
direct and diffuse solar radiation respectively. This layer is treated as a thin 
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lid covering the earth-atmosphere system in order to compute the albedo sensi- 
tivity. An assumption made here is that the reflection from the portion of the 
atmosphere above tlie layer is negligible. The altered albedo of the underlying 
system wMch is composed of the atmosphere, clouds and land is obtained from 
the multiple reflection model, Equation (4) 


a'(Mo) = a(Mo) + 


a'(Mo)[l - a(Mo)3 (1 - a) 
1 - aa'(Mo) 


(13) 


where a'( Mq ) unperturbed system albedo and a'(jUo) i® the albedo in the pres- 

ence of the aerosol layer. As there is no absorption in the model aerosol layer at 
0.55 nm (co= 1), Tc = 1 - a(jUQ) and = 1 - a. Equation (13) may be applied to 
each latitude belt and the resulting albedo integrated over each hemisphere using 
Equation (2), to yield the albedo sensitivity. Each latitude belt is identified with 
a different mean solar zenith angle, so the planar reflectivity of tlie layer, a(Mo) 
Avill be different. Also, land form and cloud cover are variable over the globe, 
so the radiative perturbation caused by the introduction of an enhanced stratos- 
pheric aerosol layer wll be latitude dependent. If instead of mean aimual aver- 
ages, seasonal effects are eonsideredj the changing solar zenith angle, cloud 
cover and land reflectivity will produce a seasonal dependence as was sliowi 
by Herman, et al. , (1976). In the next section the significance of this latitudinal 
variation will be discussed in terms of the planetary albedo sensitivity. 
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7. PLANETARY ALBEDO CHANGES 


As mentioned in Section 2, the hemispheric planetary albedo is calculated 
by integrating the albedo of each latitude belt, “p , over the hemisphere using 
Equation (2). The unperturbed zonal albedo, “p , is calculated using Equation 
(1) in which mean annual cloud cover fractions may be used. To compute the 
perturbed planetary albedo, previous authors have used several options. 
Harshvardhan and Cess (1976) assumed a system albedo of 0. 3 to calculate 
enhanced solar reflectivity; Herman, et al. , (1976) allowed for latitudinal vari- 
ation of surface albedo but assumed 50% cloud cover with a cloud top albedo of 
0. 5; and Pollack, et al. (1976) assumed a surface albedo of 0. 1 and 50% cloud 
coyer with a specified cloud deck between 3 and 6 km. In this study the mean 
annual cloud cover is not 50% over each latitude zone and the albedos of the 
cloudy and cloud-free portions of the zone are significantly different and so may 
be expected to be perturbed to different extents with the addition of a reflecting 
layer. Moreover, lower latitudes cover more area and have greater insolation 
over the year, so planetary albedo changes will be weighted in favor of low 
latitude changes. 

Apart from the above, the wide variation in planar reflectivity of the strato- 
spheric aerosols, sl(/Xq), with latitude will be a dominant effect. Table 6 gives 
the planar albedo for model aerosol layers with two scattering phase functions 
at the annual mean solar zenith angle corresponding to the different latitude belts. 
Owing to the forward bias of the phase function with g = 0. 73, the planar albedo 
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Table 6 


Albedo of Aerosol Layer 


Latitude 

Cosine of 
Mean Solar 
Zenith Aiigle 

a(Mo)/r 

g = 0. 726 

o 

11 

bfl 

N 85° 

0.257 

1.045 

1.946 

75° 

0.274 

0.941 

1.825 

65° 

0.316 

0.739 

1. 582 

55° 

0.385 

0.520 

1.299 

45° 

0.453 

0.384 

1.104 

35° 

0.513 

0.302 

0.975 

25° 

0.560 

0. 254 

0.893 

15° 

0.592 

0.227 

0,845 

N 5° 

0.609 

0.214 

0. 821 

a 

T 

0.401 

1.0 


is much less than that for die isotropic phase function. Also, there is a much 


wider variation with latitude. The effect of latitude on the albedo sensitivity may 
be illustrated by using Equation (13) to compute the change in the albedo of the 
underlying system for arbitrary unperturbed system albedos. This is presented 
in Figures 6 and 7 for a non-absorbing aerosol layer with an asymmetric phase 
function (g = 0. 73) and an isotropic phase function ( g = 0) respectively. The 
ordinate is the albedo sensitivity, here 

r 

and a' ( m) is the unperturbed system albedo. Results are presented for direct 
beam radiation at the solar zenith angles corresponding to different latitudes 
and for diffuse radiation. Two features may be noted from these figures* For 
g = 0. 73, the albedo sensitivity for latitudes equatorward of 45° is less than that 
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ALBEDO SENSITIVITY 



SURFACE REFLECTIVITY 

Figure 6. Albedo Sensitivity for Model Aerosol Layer with g = 0. 73 
for Direct-Beam Radiation at Annual Mean Solar Zenith Angles Corre- 
sponding to Latitudes 85°, 65°, 55°, 35°, 5°, and for Diffuse Radiation 
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ALBEDO SENSITIVITY 


2.0 



SURFACE REFLECTIVITY 

Figure 7. Same as Figure 6 Except with g = 0 (Isotropic Scattering) 


28 



for diffuse radiation. These latitude belts contain over 70% of the surface area 


and receive more solar insolation than regions poleward of 45°. This explains 
some of the discrepancies in planetary albedo sensitivity obtained by using dif- 
ferent averaging procedures as described later. The other notable feature is 
the almost linear relationship between the albedo sensitivity and surface reflec- 
tivity at high latitudes. This feature is present in the case of slightly absorbing 
aerosols as shown in Figure 8 in which results at two mean solar zenith angles 
corresponding to 85° N and 5°N latitude are presented for co = 0.99, 0.95 and 0.9. 


Here three methods of estimating planetary albedo changes are compared. 
For a particular set of albedos, and and cloud cover fraction C, the 
planetary albedo is computed. Now it is assumed that a miiform optically thin 
layer of aerosol covers the planet. Using tlie multiple reflection model, the 
new system albedo is 


A 


new 


^old 

= a + 

I ■ 


(14) 


Note that the diffuse or Bond albedo of the stratospheric layer, a, has been used, 
which for the present model is 0. 4r , such that the albedo sensitivity is 


^new •^old _ 0-4(1 ^old^” 
r 1 - 0.4rAo]d 


0,4(1 - AqI^j)" 


(15) 


when r«l. Tins is a hemispheric amiual averaging procedure. Another 
approach is to compute the change in reflectivity of each latitude belt with 
constant cloud cover. Tins is a zonal annual averag'e accomplished by applying 
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tlie reflection scheme to each etp and then integrating over the hemisphere using 
Equation (2). The change in planetary albedo is then 


4 f'"'" 

Ajiew “ ^old ~ ^ J ^^p,ne\v “ 0ip Q\^)Qcoii(j)d(!) 


(16) 


where, for each latitude belt 


^p,new(^) 


m + 


«p,old(l - a)[l - 
1 -a^p^old 


(17) 


with (^) the latitude. The zenith angle dependence of the planar reflectivity has 
been translated to a latitude dependence. The mean annual solar zenith angle 
has been chosen for tMs purpose. 

The zonal annual average computation may also be carried out separately 
for the cloudy and clear sl^y portions of the atmosphere. Here the multiple 
reflection model takes the form 




«c,s,oldO-a)[l-a((^))] 


1 - aa 


(18) 


c,s,old 


and 


<^p,nevv “p,old 


~ (“c.new “CjOld)^ 


new 


«s,old)n 


(19) 


The three methods of computing the albedo sensitivity have been carried 
out for the clear sky albedos and the three model cloudy albedos given in Figures 
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Table 7 

Mean Annual Albedos and Cloud Cover Fractions 


Latitude 

Cloud 

Fraction 



Model A 

Model B 

Model C 

N 85° 

0.55 

0.58 

0.67 

0.67 

0.66 

75° 

0.61 

0.54 

0.66 

0,65 

0.65 

65° 

0.64 

0.30 

0.56 

0. 53 

0.51 

55° 

0.64 

0.24 

0.52 

0.49 

0.46 

45° 

0.57 

0.20 

0. 50 

0.47 

0.43 

35° 

0.47 

0.17 

0.47 

0.44 

0.41 

25° 

0.41 

0.17 

0.46 

0.43 

0.40 

15° 

0.44 . 

0.13 

0.44 

0.41 

0.38 

N 5° 

0.51 

0.12 

0.44 

0.40 

0. 36 

S 5° 

0.50 

0.12 

0. 44 

0.40 

0.36 

15° 

0.47 

0.13 

0.44 

0.41 

0. 37 

25° 

0.47 

0.14 

0.45 

0. 42 

0. 38 

35° 

0.54 

0.14 

0.47 

0.43 

0, 39 

45° 

0.65 

0.15 

0.48 

0.45 

0.41 

55° 

0.79 

0.18 

0.50 

0.47 

0.43 

65° 

0.77 

0.26 

0.54 

0.52 

0.49 

75° 

0.56 

0.60 

0,69 

0, 69 

0.66 

S 85° 

0.47 

0.61 

0.70 

0.69 

0.69 
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Table 8 

Albedo Sensitivity 


Models 

Planetary 

Albedo 

Albedo Sensitivity a A/at 

Hemispheric 
Annual Average 

Zonal Annual 
Average (Cloudy 
& Clear 
Combined) 

Zonal Annual 
Average (Cloudy 
& Clear 
Separate ) 

N.H. 

S. H. 

N.H. 

S. H. 

N. H. 

S. II. 

N. H. 

S. H. 

Model A 

0. 323 

0. 327 

0. 184 

0. 182 

0. 123 

Q. 123 

0. 132 

0. 133 

Model B 

0.308 

0.309 

0. 192 

0. 191 

0. 130 

0. 132 

0. 137 

0. 139 

Model C 

0.291 

0.289 

0. 202 

0, 203 

0. 138 

0. 142 

0. 143 

0. 1^7 

Cess 

0. 306 

0. 301 

0. 193 

0. 196 

0. 132 

0. 136 



0. 138 

0. 140 




2 and 4 respectively. The mean annual albedos and cloud cover fractions from 
London (1957) and van Loon (1972) for the N and S Hemispheres respectively 
are given in Table 7. The hemispheric albedo sensitivity calculated using the 
three methods are listed in Table 8 along with the miperturbed planetary albedo 
for the three different cloud models and also for albedo values deduced by Cess 
(1976) from satellite derived results in Ellis and Vender Haar (1976). 

The results of the albedo sensitivity computation are significant in two ways. 
First, there is a substantial reduction in the sensitivity when a zonal average is 
taken instead of a hemispheric average. The mean value of tliis reduction is 
31% for the four cases considered. Tliis lo^ver value is a more appropriate 
measure of the mean annual planetary albedo sensitivity to a muform strato- 
spheric aerosol perturbation. The explanation for this reduction may be fomid in 
Table 6 and Figure 6. The planar reflectivity of the aerosol layer is a very 
strong function of the zenith angle and the use of hemispheric averaging is equiv- 
alent to assmning a diffuse reflectivity wliich is greater than the planar reflec- 
tivity for mean solar zenith angles corresponding to latitudes equatorward of 
45°, The increased reflectivity of the layer in polar regions is compensated by 
the effect of lower latitudes wluch have more area and receive more solar in- 
solation. If the aerosol scattering had been isotropic this effect would not have 
been as pronounced. Again from Table 6 and Figure 7, it is evident that tlie 
planar reflectivity for isotropic scattering is nearly miiform till 35° latitude and 
is very close in value to the diffuse reflectivity. An example is a computation 
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as presented in Table 8 for isotropic scattering. Using the parameters of Model 
B for the Northern Hemisphere, the albedo sensitivity is 0,480 with hemispheric 
a,veraging and 0,465 with zonal averaging (0,482 when cloud and cloud-free 
albedos are taken separately). 

The other feature of the results is the insignificant change in albedo sensi- 
tivity when cloudy and clear sky portions of each latitude zone are considered 
separately. This may seem surprising because the albedo sensitivity is a strong 
function of the reflectivity of the underlying surface as is shown in Figure 6, 
However, the relationship is monotonic and nearly linear for the higher latitudes; 
the same is true for lower latitudes when the underlying system albedo is less 
than 0. 5. The effective albedo of a partly cloudy area lies between the cloudy 
and clear-sky albedos, the exact value depending on the cloud cover fraction, 

C, Equation (1). Owing to the near linearity of the albedo sensitivity relation- 
ship, the value obtained by considering cloudy and clear-sky portions separately 
and considering an effective zonal albedo are very nearly the same. The sig- 
nificance of this result is that it is not essential to know the exact cloud cover 
fraction and cloudy and clear-sky albedos to compute the albedo sensitivity in 
this case. This argument holds for slightly absorbing aerosols as well, as show 
shown in Figure 8, Satellite measurements of effective zonal albedos are avail- 
able by month in 10° latitude belts (Ellis and Vender Haar, 1976) and present 
ERB studies are' geared to obtaining this information. Therefore it is possible 
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Figure 8. Same as Figure 6 Except for only 85° and 5° Latitude 
with CO = 1,0, 0.99, 0.95 and 0.9 
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to estimate the albedo sensitivity of different latitude zones over the aimual cjmle 


with current data. 

8. MONTHLY ALBEDO SENSITIVITY 

The procedure for computing the monthly albedo sensitivity is similar to 
the annual case. The imperturbed zonal albedo, aoid , is taken to be the mean 
monthly values quoted in Ellis and Vender Haar (1876). By the discussion in the 
previous section, we need not take into accomit the cloudy and cloud-free portions 
of the zone separately. The reflectivity of the aerosol layer, a(pQ), is a function 
of the solar zenith angle, the variation of which must be included in the compu- 
tation. To simplify matters, the 21st day of each month has been chosen to 
compute the reflectivity wMch has been averaged over the zenith angles made 
by the path of the sun on that day for the mid-point of each latitude belt. The 
albedo perturbation caused by the presence of a non-absorbing stratospheric 
aerosol layer of incremental optical triickness 5 Ar , vlll then be {^new " ^'old ) 
where a^ew is given by Equation (17). With the assumption that the optical depth 
at 0. 55 ^m (denoted by Tyis) is representative of the aerosol layer, the albedo 
sensitivity, Aa /Ar , obtained from Equation (17) has been plotted in Figure 9. 

The isopleths in Figures 9 and 10 are in units of change of albedo per unit 
change in optical depth or alternately, the change in albedo in percent of incident 
solar energy for an optical depth perturbation of Ar = 0. 01. The dominant 
feature is seen to be the zenith angle dependence translated here in a monthly 
and latitudinal dependence. Over high latitudes the effect of surface reflectivity 
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SOUTH LATITUDE north 










may also be noticed, such as the lower values of Act/Ar in late winter and early 
spring when there is snow cover on the ground. This phenomenon has been 
discussed in Section 7 with the help of Figure 6. Figure 10 is a plot of tlie solar 
averaged albedo sensitivity which is obtained by repeating tlie computations 
carried out at 0. 55 juni at a number of wavelengths and taldng a weighted mean. 
The aerosol considered here is slightly absorbing in the solar infrared and 
Equation (17) may be rephrased in terms of the transmission as 

«new = a + (20) 

where T is planar and T is spherical transmission. 

These figures may be compared with Figures 5 and 11, respectively, 
of Herman, et al. (1976) who have plotted the change in albedo in percent 
for a stratospheric aerosol perturbation optical depth of 0.03 at 0.5 (xm. 
Their results are for the Northern Hemisphere for a prescribed surface albedo 
variation which assigns values of the albedo according to tliree classes of sur- 
face (bare ground, open water and ice or snow). Clouds are assumed to cover 
50% of the area and are simulated by results corresponding to a ground albedo 
of 0. 5. Their results for non-absorbing aerosols show veiy clearly the high- 
latitude regions of minimum albedo sensitivity in spring and maximum albedo 
sensitivity in tiie fall. Their results for 75% H7SO4 averaged over solar wave- 
lengths correspond quite closely to Figure 10 presented here except for Mgh 


39 


LATITUDE 



Figure 11. Change in System Albedo in 1 ercent Nsitli Addition 
of Aerosol Layer of Optical Depth 0.03 at 0. 5/im: Present 
Model (Above) and Herman, et al. (1976) (Below) 
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latitude regions in the spring and fall. The comparison is facilitated by replotting 
Figure 10 for the Northern Hemisphere for an opticEil depth of 0. 03 at 0. 5 [xm. 

This is shown in Figure 11 in wliich Figure 11 of Herman, et al. , (1976) has been 
reproduced below for comparison. (It should be noted that the scattering phase 
function used in their work is slightly different from the one used here. ) One 
reason for the discrepancy is the low values of albedo for high latitude regions 
quoted by Ellis and Vonder Haar (1976), who recognize that their results are 
quite imcertain in the 70°- 90° latitude range. 

9. RADIATION BALANCE AND CLIMATE MODELING 

The annual and monthly zonal albedo sensitivity results have important im- 
plications in energy -balance climate modeling. The sensitivity of these models 
is usually related to perturbations in the magnitude of the solar constant. For 
example, both Sellers (1969) and Budyko (1969) refer to an increase in volcanic 
dust loading and a decrease in the solar constant interchangeably. However, 
even a uniform dust layer of global extent will not lead to a uniform reduction 
in the solar energy received by a colunm of the earth-atmosphere system below 
the stratospheric layer. The computations presented here have shown how tliis 
reduction is distributed in latitude and with season. The basic equation of zonal 
energy balance on an amiual mean basis is (Budyko, 1969) 

Q(0jTl - a(0)] - 1(0) = D(0) (21) 
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where Q(0) is the solar insolation at latitude 0, a{(p) the albedo, 1(0) is the out- 
going terrestrial infrared radiation and D(0) is the divergence of heat. The 
presence of a stratospheric aerosol layer modifies a through increased reflec- 
tivity and I through an increased greenhouse effect acting in the opposite direc- 
tion. The reduction in solar energy absorbed, QAct, caused by a uniform 
stratospheric layer of Tyis = 0. 1 is plotted in Figure 12 which follows from 
Figure 10 and the distribution of Q(0) for each month. 

The net loss in a column of the earth- atmosphere system, however, is 
QAa- AI, which is the change in the radiation balance and it is necessary to 
compute the infrared greenhouse sensitivity of the aerosol layer. This is done 
using an atmospheric radiation model similar to the one by Harshvardhan and 
Cess (1978) in which an emissivity formulation has been used to calculate the 
infrared flux. Water vapor and carbon dioxide are considered to be the gaseous 
absorbers in the atmosphere and a 4km thick layer of the model aerosol layer 
is added to determine the enhanced greenhouse effect. The computations have 
been carried out for the model atmospheres given by McClatchey, et al. , (1971) 
for clear skies and a cloudy sky with cloud top at 6. 5 km. The aerosol layer 
was positioned to correspond to the observed level of peak concentration at dif- 
ferent latitudes (Rosen, et al. , 1975). Results are presented in Table 9 in which 
the infrared greenhouse sensitivity is expressed in (W/m^)/r^j^ and is seen to be 
a maximum in the tropics and a minimum for the subarctic winter model. 
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Figure 12. Reduction in Mean Monthly Solar Energy Absorbed 
in W/m with Addition of Aerosol Layer of r . = 0, 1 
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Table 9 

Infrared Greenhouse Sensiti'^dly 


Model 

Infrared Greenhouse Sensitivity (W/m^)/Tvji, 

Clear Sky 

Cloudy Sky 

Tropical 

14.8 

7.5 

Midlatitude 

13.8 

6.9 

Summer 

Midlatitude 

8.6 

3.6 

Winter 


Subarctic 

10.5 

4.0 

Simimer 

Subarctic 

5.6 

1.8 

Winter 


An increase in the aerosol optical thiclmess of Tyjg - 0. 1, for example, would 
result in an increase in radiation trapped in the earth-atmosphere column under 
clear sky conditions of 1. 5 W/m~ in the tropics and 0. 6 W/m“ in the subarctic 
winter. Cloudiness may be taken into account by a weighting of the clear and 
cloudy sky results. This greenliouse moderation, when combined with solar 
depletion results in the distribution of net energy loss, QAq' - AI as shown in 
Figure 13 for Tyjs = 0.1. Rather than use climatological cloudiness, the average 
of the clear and cloudy sky results has been used. Also, the greenhouse effect 
in spring and fall have been taken as intermediate between winter and summer. 
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Figure 13. Reduction in the Zonal Radiation Balance in W/m^ 
with Addition of Aerosol Layer of - 0. i 
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One observation from Figure 13 is that the presence of a stratospheric aero- 
sol laj^er composed of 75 % H 2 SO 4 results in a loss of radiation absorbed within 
the earth- atmosphere system at all latitudes and seasons except the polar winter. 
The relationsMp between tliis net decrease in energy absorbed and mean atmos- 
pheric or surface temperature may be established through climate models. The 
maximum perturbation in the radiation balance occurs at high latitudes in the 
spring and fall whereas polar winters show a net gain in the radiation balance due 
to the increased greenhouse effect. Latittides equatorward of 50° show a nearly 
uniform reduction in the radiation balajace. 

The magnitude of the change in the net radiation is only a small percentage 
of the range of typical values of the radiation balance.?- Monthly radiation balance 
over all latitudes have been presented by Ellis and Vender Haar (1976) and values 
range from -175 W/m- in the south polar winters to +110 W/ni- in the southern 
hemisphere tropics in summer. The uncertainty in mean montlily net radiation 
measurements ranges from 3 to 14W/m^ such that the computed values of the 
change in the radiation balance caused by the aerosol layer of ryjg = 0. 1 are less 
than the uncertainty in the satellite derived data. It is unlikely that there \vill be 
a greater perturbation to the net radiation for extended periods of time following 
volcanic eruptions. Following the Agung eruption, Tyjs 0. 2 in the Southern 
Hemisphere for a year (Volz, 1970) and even eruptions of greater magnitude 
will not alter the net radiation appreciably. 
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The perturbation in the radiation balance can also be expressed as a change 
in the effective black-body radiative temperature of die earth-atmosphere sys- 
tem. Tills change is given by 


AT 


QAa - AI 
4uT3 


( 21 ) 


where a is the Stefan- Boltzmaim constant, T the effective black-body radiative 
temperature of the earth- atmosphere system and (QAa - A I) the change in the 
radiation balance. Figure 14 is a plot of AT in °K for = 0. 1. The major 
features detected in the plot of radiation balance perturbation appear enhanced 
in tills representation. It can be seen that the major radiative effect of the 
stratospheric aergsol layer occurs in the spring and fall when there is a pro- 
nounced change in the equator to pole radiative energy gradient. As the results 
pertain to the entire earth-atmosphere column, it is not possible to estimate the 
change in the diabatic heating at various levels of the a^tniosphere or estimate 
surface temperature changes. However, it axipears from tliis mialysis that a 
uniform layer of stratospheric aerosols would have only a small effect on the 
long term radiative regime equr.torwards of 50°~60° if the change in turbidity 
corresponded to that caused by the Agimg eruption. The localized radiative 
effect in the vicinity of the aerosol layer could still be considerable as discussed 
in Harshvardhan and Cess (1976). 
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LATITUDE 



Figure 14. Change in the Effective Black- Body Radiative 
Temperature of the Earth -Atmosphere System in K with 
Addition of Aerosol Layer of = 0. 1 
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10. SUMMARY 


The present study has examined the role of stratospheric aerosols in altering 
the radiation balance of the earth-atmosphere system. It has been confirmed that 
for an aerosol layer of 75% H 2 SO 4 the dominant influence is an increased re- 
flection of solar energi^ all over the globe except for tlie polar winters when tliere 
is no solar insolation. An important conclusion resulting from the albedo sensi- 
tivity computations is the necessity of including solar zenith angle effects in 
considering albedo enhancement. It was also pointed out that die separation of 
areas into cloud-covered and cloud-free is unnecessary in the calculation of the 
radiative effects of the reflecting layer. This enables the use of mean monthly 
planetary albedos obtained from satellite measurements to calculate the change 
in tlie radiation balance of each zonal belt caused by the presence of a strato- 
spheric layer. The xierturbation in the radiative balance may be used to simu- 
late the radiative effect of a stratospheric aerosol layer in climate studies. 
Results indicate that the perturbation is strongest at liigh latitudes, particularly 
in the polar region, in the spring and fall. 



REFERENCES 


Braslau, N. and J. V. Dave, 1973 ; Effect of aerosols on the transfer of solar 
energy through realistic model atmospheres. Part II: partly absorbing 
aerosols. J. Appl. Meteor., Vol. 12, pp. 616-619. 

Budyko, M. I. , 1969: The effect of solar radiation variations on the climate of 
the earth. Tellus , Vol. 21, pp. 611-619. 

, 1974: Climate and Life. Academic Press, New York, 508 pp. 

Cadle, R. D. and G. W. Grams, 1975: Stratospheric aerosol particles and their 
optical properties. Rev. Geophys. Space Phys., Vol. 13, pp. 475-501. 

Cess, R. D. , 1976: Climate change: an appraisal of atmospheric feedback 
mechanisms employing zonal climatology. J. Atmos. Sci. , Vol. 33, pp. 
1831-1843. 

Coakley, J. A. and P. Chylek, 1975: The two-stream approximation in radiative 
transfer; including the angle of the incident radiation. J. Atmos. Sci. , 

Vol. 32, pp. 409-418. 

Curran, R. J. , R. Wexler and M. L. Nack, 1978: Albedo climatology analysis 
and the determination of fractional cloud cover. NASA Tech. Memo. 

79576. 


50 



Dave, J. V. and N. Braslau, 1974: Effect of cloudiness on the transfer of solar 
energy through realisitc model atmospheres. Kept. RC 4869, IBM Thomas 
J. Watson Research Center, Yorktown Heights, N. Y. 

and , 1975: Effect of cloudiness on the transfer of solar 

energy through realistic model atmospheres. J. Appl. Meteor., Vol. 14, 
pp. 388-395. 

Dyer, A. J. and B. B. Hicks, 19'68: Global spread of volcanic dust from the 
Bali eruption of 1963. Quart. J. Roy. Meteor. Soc., Vol. 98, pp. 545-554. 

Ellis, J. and T. H. Vonder Haar, 1976: Zonal average earth radiation budget 
measurements from satellites for climate studies. Atmos. Sci. Pap. 240, 
Colorado State University, Fort Collins, Co. 

Hansen, J. , 1969: Exact and approximate solutions for multiple scattering by 
cloudy and hazy planetary atmospheres. J. Atmos. Sci. , Vol. 26, pp. 
478-487. 

Harshvardlian and R. D. Cess, 1976: Stratospheric aerosols : effect upon 

atmospheric temperature and global climate. Tellus , Vol. 28, pp* 1-10. 

Herman, B. M. , S. R. Browning, and R. Rabinoff, 1976: The change in earth- 
atmosphere albedo and radiational equilibrium temperatures due to strato- 
spheric pollution. J, Appl. Meteor., Vol. 15, pp. 1057-1067. 


51 




Hofmann, D. J. , J. M. Rosen, T. J. Pepin, and R. G, Pinnick, 1975; Strato- 
spheric aerosol measurements I: time variations at Northern midlatitudes. 
J. Atmos. Sci. . Vol. 32, pp. 1446-145G. 

Jacobomtz, H. , W. L. Smith, il. B. Howell, F. W. Nagle, and J. R. Hickey, 
1978: The first eighteen months of planetary radiation budget measure- 
ments from the Nimbus 6 ERB experiment. Preprints Third Conf. 
Atmospheric Radiation, Davis, Ca. , Amer. Meteor. Soc. , pp. 164-166. 

Kondratyev, K. Ya. (Ed.), 1975: Radiation characteristics of the Atmosphere 
and the Earth’s surface. NASA TTF-678, Amerind Publishing Co. , New 
Delhi, 580 pp. 

London, J. , 1957; A study of the atmospheric heat balance. Report, Contract 
AF19(122)-165, College of Engineering, New York University. 

Luther, F. M. , 1976: Relative influence of stratospheric aerosols on solar and 
longwave radiative fliixes for a tropical atmosphere. J. Appl. Meteor. , 
Vol. 15, pp. 951-955. 

McClatchey, R. A., R. W. Fenn, J. E. A. Selby, F. E. Volz, and J. S. 

Garing, 1971: Optical propei'ties of the atmosphere (rev.), Environ. Res. 
Pap. 54, Air Force Cambridge Res. Lab., Hanscom Air Force Base, 
Bedford, Mass. 


52 



Palmer, K. F. and D. Williams, 1975: Optical constants of sulfuric acid; 
application to the clouds of Venus ? Appl. Optics, VoL 14, pp. 208-219. 

Pinnick, R. G. , J. M. Rosen, and D. J. Hofmann, 1976: Stratospheric aerosol 
measurements III: optical model calculations. J. Atmos. S ci. , Vol. , 33, 
pp. 304-314. 

Pollack, J. B. , O. B. Toon, C. Sagan, A, Summers, B. Baldmn, and W. 
van Camp, 1976: Volcanic explosions and climatic change: a theoretical 
assessment. J. Geophys. Res., Vol. 81, pp. 1071-1083. 

Rosen, J. M. , 1971: The boiling point of stratospheric aerosols. J. Appl. 
Meteor. , Vol. 10, pp. 1044-1046. 

, D. J. Hofmann, and J. Laby, 1975: Stratospheric aerosol measure- 
ments II: the worldmde distribution. J. Atmos. Sci. , Vol. 32, pp. 1457- 
1462. 

Russell, P, B. and R. D. Hake, Jr,, 1977: The post- Fuego stratospheric aero- 
sol: lidar measurements, with radiative and thermal implications. J. 
Atmos. Sci. , Vol. 34, pp. 163-177. 

Sellers, W. D. , 1969: A global climatic model based on the energy balance of 
the earth-atmosphere system. J. Appl. Meteor. , Vol. 8, pp. 392-400. 


53 



q ■:-;.rL...:zJ J.... -. J '^fL i. 

Shettle, E. P. and B,. W. Fenn, 1976: Models of the atmospheric aerosols and 
their optical properties. AGAE-D Conference Proceedings No. 183 on Optical 
Propagation in the Atmosphere. 

van Loon, H. , 1972: Cloudiness and precipitation. Meteor. Monogr. , No, 35, 

pp. 101-111. 

Volz, F. E. , 970: Atmospheric turbidity after the Agung eruption of 1963 and 

size distribution of the volcanic aerosol. J. Geophys. Res. , Vol, 75, 
pp, 5185-5193, 

Vender Haar, T. H. and J. Ellis, 1975: Albedo of the cloud-free earth- 
atmosphere system. Preprints Second Conf. Atmospheric Radiation, 
Arlington, Va. , Amer, Meteor. Soc., pp. 107-110. 

Wiscombe, W. J. , 1975: Solar radiation calculations for arctic summer stratus 
conditions. Climate of the Arctic, G. Weller and S. A. Bowling, Eds. , 
Geophysical Institute, University of Alaska, pp. 245-254. 

and G. W. Grams, 1976: The backscattered fraction in two-stream 

approximations, Atmos. Sci., Vol. 33, pp. 2440-2451. 


54 


